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Abstract 



^ . In the complex plane, the frequency response of a univariate polynomial is the 

\ set of values taken by the polynomial when evaluated along the imaginary axis. This 

is an algebraic curve partitioning the plane into several connected components. In 
this note it is shown that the component including the origin is exactly representable 
by a linear matrix inequality if and only if the polynomial is stable, in the sense 
O ' tliat all its roots have negative real parts. 
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1 Introduction 



Let G M[s] be a polynomial of the complex indeterminate s G C. We say that is 
stable if all its roots lie in the open left half-plane. Define the frequency response 

V = {p{juj) : ooeR} 

as the set of values taken by the polynomial when evaluated along the stability boundary, 
namely the imaginary axis. The frequency response plays a key role when deriving results 
of robust control theory such as Kharitonov's theorem [21 [Sj |5]. 

In [H] it was observed that the frequency response of a stable polynomial features inter- 
esting convexity properties, see also [3^, Chapter 18]. More specifically, given a polynomial 
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p{s), an arc is defined as a subset of tlie frequency response for a given range of the inde- 
terminate, i.e. {p{juj) : uj & Q} with Q a subset of the real hne M. A proper arc is an arc 
that does not pass through the origin and such that the net change in the argument of 
p{joj) does not exceed vr as a; increases over VL. The Arc Convexity Theorem of [9] states 
that all proper arcs of the value set of a stable polynomial are convex. Alternative proofs 
can be found in [TT] and [S]. 

The frequency response dV is a curve that partitions the complex plane into several 
connected components. We denote by V the connected component including the origin. 
It is called the inner frequency response set in [9] . In the case of a stable polynomial, the 
boundary of V therefore consists of a finite union of proper arcs. Theorem 4.1 in [9] uses 
the Arc Convexity Theorem to establish convexity of V. 

In this note we provide a more accurate description of the geometry of V and an alternative 
proof of its convexity. We derive an explicit representation of this set as a two-dimensional 
linear matrix inequality (LMI). Instrumental to this derivation are standard results from 
real algebraic geometry [H IH O [H] and a recent characterization of two-dimensional 
convex polynomial level sets obtained in jlU]. 



2 Algebraic description of the frequency response 



The frequency response of polynomial p{s) = Pq + PiS + ■ ■ ■ + PnS"" can be expressed as 
the parametric curve 

dV = {pijuj) = x{iu) +]y{to) : G M} (l) 

where 



x{uj) 

qy{u 



(2) 



and 

q^iu) = po- P2UJ^ + Pauj^ H 

qy{uj) = piuj - p:iUJ^ + p^u^ ^ 

q^{uj) = 1 

are polynomials of the real indeterminate a; G M. 

Curve dV is rationally (here polynomially) parametrized, so it is an algebraic plane curve 
of genus zero In control theory terminology, dV is sometimes called the Mikhailov plot 
of polynomial p{s) or the Nyquist plot of the rational (here polynomial) transfer function 
p{s), see [3] or [S]. 

Equations ([lH2]) provide a parametric description of curve dV. With the help of elimina- 
tion theory and resultants, we can derive an implicit description 

dV = ix+w : f(x,v) = 01 (31 



where f{x,y) is an irreducible bivariate polynomial, see [HI Section 3.3]. 



Lemma 1 Given two univariate polynomials g{uj) = go + QiUJ + ■ ■ ■ + Qn^^"' £ M.[uj] and 
h{uj) = hQ + hiUJ + - ■ ■ + hnUj" G W[uj], there exists a unique (up to sign) irreducible polyno- 
mial b{g, h) G ffi[5'0) fl'i) • • • 5 gn, ho, hi, ... , called the resultant which vanishes whenever 
g{uj) and h{uj) have a common zero. 



To address the implicitization problem, we make use of a particular resultant, the Bezoutian, 
see e.g. [Til Theorem 4.1] or [7^ Section 5.1.2]. Given two univariate polynomials 
g{uj),h{uj) of the same degree n as in Lemma [H build the following bivariate polyno- 
mial 

g{u^)h{v) - g{v)h{u) ^ ^ ^ 

UJ — V ^ ^ 

k=0 1=0 

called the Bezoutian of g and h, and the corresponding symmetric matrix B{g, h) of size 
n X n with entries bki bilinear in coefficients of g and h. As shown e.g. in [lU Theorem 
4.3] or [TJ Section 5.1.2], the determinant of the Bezoutian matrix is the resultant: 



Lemma 2 det B{g,h) =b{g,h). 

Now we can use the Bezoutian to derive the implicit equation ([3]) of curve dV from the 
explicit equations ([I]l2]). 

Lemma 3 Given polynomials qx,%, Iz of equations (12113', polynomial of equation Iji^ 
is given by f{x,y) = det F{x,y) where 

F{x, y) = ±{B{qx, qy) - xB{qy, q,) - yB{qx, q^)) (4) 

is a symmetric pencil, i.e. a polynomial matrix linear in x, y, of size n. 



Proof: Rewrite the system of equations ^ as 

g{uj) = qx{u) - xqz{uj) = 
h{uj) = qy{uj) - yqziuj) = 

and use the Bezoutian resultant of Lemma 12 to eliminate indeterminate u and obtain 
conditions for a point {x,y) to belong to the curve. The Bezoutian matrix is B{g,h) = 
B{q^ - xq^, qy - yq^) = B{q^, q^) - xB{qy, q^) - yB{q^, q^). Linearity in x, y follows from 
bilinearity of the Bezoutian. Finally, note that the sign of F{x,y) = ±B{g, h) affects the 
sign of /, but not the implicit description f{x,y) = O.D 

Lemma [3] provides the implicit equation of curve dV in symmetric linear determinantal 
form. 



3 Convexity properties of the inner frequency re- 
sponse set 



Curve dV partitions the complex plane into several connected regions. We are interested 
in the connected region containing the origin, denoted by V. In order to study the 
geometry of this region, we need the following result. 

Lemma 4 The sign of pencil F{x,y) in ^ can be chosen such that F{0,0) = B{qy,qz) 
is positive definite if and only if p is a stable polynomial. 

Proof: The signature of the Bezoutian matrix B{qx,qy), (the number of positive eigen- 
values minus the number of negative eigenvalues) is equal to the Cauchy index of the 
rational function g^, (a;)/ 5^,(0;) (the number of jumps from —00 to +00 minus the number 
of jumps from +00 to —00), see [H Section 9.1.2]. The Cauchy index is maximum (resp. 
minimum) when B{qx,qz) is positive (resp. negative) definite. This occurs if and only if 
polynomials qx and qy satisfy the root interlacing condition, i.e. they must have only real 
roots and between two roots of qxi^^) there is only one root of qy{uj) and vice- versa. Since 
qx{uj) = Rep{juj) and qy^uj) = lmp{ju), this is equivalent to stabihty of p in virtue of the 
Hermite-Biehler theorem, see pi Section 8.1] or ^ Section 1.3].n 

The main result of this note can now be stated. 

Theorem 1 The connected component including the origin and delimited by the frequency 
response of polynomial p can be described by a linear matrix inequality (LMI) 

V = {x+jy : F{x,y)hO} 

if and only if p is stable. In the above description, F{x,y) is given by and >z means 
positive semidefinite. 

Proof: First we prove that p stable implies LMI representability of V. This set is the 
closure of the connected component of the polynomial level set {x,y : f{x,y) > 0} that 
contains the origin, an algebraic interior in the terminology of [TO]. Polynomial f{x,y) is 
called the defining polynomial. By continuity, the boundary of V consists of those points 
x + jy for which F{x, y) drops rank while staying positive semidefinite. Note that in 
general this boundary is only a subset of the curve dV. 

To prove the converse, namely that LMI representability of V implies stability of p, we use 
a result of [10] stating that a two-dimensional algebraic interior containing the origin has 
an LMI representation if and only if it is rigidly convex. Geometrically, this means that 
a generic line passing through the origin must intersect the algebraic curve f{x,y) = 
a number of times equal to the degree of f{x,y). Rigid convexity implies that x{u) 
and y{uj), the respective real and imaginary parts of polynomial p{juj), satisfy the root 
interlacing property, and this implies stability of p by the Hermite-Biehler Theorem used 
alreadv in the nroof of Lemma HI □ 



Figure 1: Frequency response of a stable third degree polynomial. The shaded region is 
the convex component containing the origin. 



4 Examples 

4.1 Stable third degree polynomial 

Let p{s) = s'^ + + 4s + 1. Then qx{uj) = 1 — a;^ and qy{uj) = — u^ in parametrization 
([2]). With the help of the Control System Toolbox for Matlab, a visual representation of 
curve dV can be obtained as follows: 

» p = [1 1 4 1] ; 7„ polynomial in Matlab format 

» nyquist (tf (p, 1) ) % frequency response 

» axis ([-5 2 -7 7]) % zoom around the origin 

see Figure [1] The Bezoutian matrices of Lemma [3] can be computed with the following 
Maple 10 instructions: 

> with (Linear Algebra) : 



> qx:=l-w"2:qy:=4*w-w"3:qz:=l: 

> Bxy :=BezoutMatrix(qx,qy ,w,method=symmetric) ; 

[-1 1] 

[ ] 

Bxy := [ -3 0] 

[ ] 
[1 0-4] 

> Byz:=BezoutMatrix(qy,qz,w,method=symmetric) ; 

[0 1] 

[ ] 

Byz := [ 1 0] 

[ ] 
[1 0-4] 

> Bxz:=subs(e=0,BezoutMatrix(qx,qz+e*w"3,w,method=syimnetric)) ; 

[ 0] 

[ ] 

Bxz := [ -1] 

[ ] 

[0-1 0] 



Note the use of the subs instruction to ensure that the last Bezoutian matrix has appro- 
priate dimension 3. Matrix B{qx,qz) is negative definite, so a sign change is required to 
build 

1 -1+x 

F{x, y) = 3 + a; —y 
— 1 + X —y A — 4x 

and we obtain the following determinantal polynomial: 



> F : =- (Bxy-x*Byz-y*Bxz) ; 

> f : =Determincint (F) ; 

2 2 c 
f:=9-3x-5x - y - x 



describing algebraic plane curve dV implicitly. The curve can be studied with the 
algcurves package of Maple: 



> with(algcurves) : 

> genus (f,x,y) ; 



> plot_real_curve(f ,x,y) ; 
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X 

Figure 2: Frequency response of a stable eighth degree polynomial. The shaded region is 
rigidly convex. 

4.2 Stable eighth degree polynomial 

A more complicated example is the eighth degree stable polynomial p{s) = 336 -|- 198s + 
496s^ + 117s^ -|- 183s'* -|- 20s^ -|- 24s^ + s^ -|- s^ whose frequency response is represented on 
Figure [2l The rigidly convex region around the origin has the LMI description 
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4.3 Unstable fourth degree polynomial 



This example is taken from [TT]. Let p{s) = s'^ — s — 1. The imphcit equation of dV is 



f{x,y) = det 
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The component V including the origin is convex, see Figure [31 but it is not rigidly convex 
since a generic line passing through the origin cuts the quartic dV only twice. Hence 
this region does not admit an LMI representation, and by Theorem [H polynomial p{s) is 
unstable. 

Note however that V can be represented as the projection of an LMI set: 



V = {x + jy : 3z : 



l+x z 


^0, 


z y 


z 1 







^0} 



by introducing a lifting variable 2;, but such constructions are out of the scope of this 
note. 



5 Conclusion 



Convexity of the connected component containing the origin and delimited by the fre- 
quency response of a stable polynomial was already established in [H]. In this note we 
give an alternative proof of this result based on Bezoutians and we give a more accurate 
characterization of the geometry of this region. Namely, the region is rigidly convex in 
the sense of [10] , a property which is stronger than convexity, and which is equivalent to 
the existence of an LMI representation of the set. 

In the terminology of convex analysis, the polynomial f{x,y) defining implicitly the fre- 
quency response in is hyperbolic with respect to the origin. Equivalently, f{x,y) can 
be expressed as the determinant of a symmetric pencil which is positive definite at the 
origin. See [13] for a tutorial on hyperbolic polynomials and [12] for connections with the 
results of [TO]. This note therefore unveils a link between polynomial hyperbolicity and 
stability. 
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Figure 3: Frequency response of an unstable fourth degree polynomial. The shaded region 
is convex, but not rigidly convex. 
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